Impurity-induced Local Density of States in a D-wave Superconductor 

Carrying a Supercurrent 
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The local density of states (LDOS) and its Fourier component induced by a unitary impurity in 
a supercurrent-carrying d-wave superconductor are investigated. Both of these quantities possess a 
reflection symmetry about the line passing through the impurity site and along the supercurrent if it 
is applied along the antinodal or nodal direction. With increasing supercurrent, both the coherence 
and resonant peaks in the LDOS are suppressed and slightly broadened. Under a supercurrent along 
the antinodal direction, the coherence peaks split into double peaks. The modulation wavevectors 
associated with elastic scatterings of quasiparticles by the defect from one constant-energy piece of 
the Fermi surface to another are displayed as bright or dark spots in the Fourier space of the LDOS 
image, and they may be suppressed or enhanced, and shifted depending on the applied current and 
the bias voltage. 



PACS numbers: 74.25.-q, 74.20.-z, 74.62. Dh 

The understanding of the local physics in cuprate or 
high temperature superconductors (HTS) is one of the 
most challenging problems in condensed matter physics 
today. Different from the conventional s-wave supercon- 
ductors, the HTS have very complex phase diagrams de- 
pending on doping and chemical composition, ft is also 
well established that the superconducting order parame- 
ter in the cuprates has predominantly d-wave symmetry 
[1]. The zero bias conductance peak (ZBCP) in the tun- 
neling spectroscopy of a normal metal-cuprate supercon- 
ductor junction with non-(nOm) contact provides one of 
the direct evidences for this symmetry [2] . Due to the d- 
wave nature of the order parameter, impurities inserted 
into cuprates can serve as an important tool to explore 
the physics of HTS. Theoretical calculations of the lo- 
cal density of states (LDOS) predicted that a strong po- 
tential scatterer could induce a resonance peak near the 
Fermi level at sites near the impurity [3,4]. This resonant 
peak near zero-bias voltage was observed at and near the 
sites of Zn impurities in Bi 2 Sr 2 CaCu 2 08 +( 5 by scanning 
tunneling microscopy (STM)[5]. In addition, an interfer- 
ence pattern with four-fold symmetry was also detected 
in the STM image [5]. 

When a superconductor carries a supercurrent (J s ), 
Cooper pairs with finite momentum appear in the sys- 
tem. This would drastically affect the electronic struc- 
ture of the superconductor including the elementary exci- 
tation spectrum, the order parameter symmetry, and the 
tunneling spectroscopy [6,7,8]. With increasing supercur- 
rent velocity, the superconducting order parameter can 
be depressed. Meanwhile, the supercurrent density first 
increases monotonously and then arrives at a maximum 
value, which is called the critical current density. Beyond 
that, superconductivity becomes unstable and collapses 
to the normal state. So the supercurrent in the stable 
regime can be also used as a probe to further understand 



the quasiparticle excitations in HTS. Moreover, a better 
understanding of the property of a superconductor under 
an applied J s may have the potential for device applica- 
tions. 

In Ref. [7] , we have studied the tunneling conductance 
characteristics between a normal metal and a d-wave su- 
perconductor (dSC) carrying a supercurrent parallel to 
the interface of the junction. It was shown that for suffi- 
ciently large applied current, the midgap-surface-state- 
induced ZBCP splits into two peaks in the tunneling 
regime. So far there exist no experimental measurements 
which could be used to compare with our theoretical pre- 
dictions. The closest tunneling experiment to the idea in 
Ref. [7] was done on YBCO under a spin injected cur- 
rent [9], it would be interesting to see that the similar 
experiment will be performed on an HTS sample car- 
rying a supercurrent in the near future. As a natural 
extension of our previous work [7], here we examine the 
LDOS induced by a strong defect which replaces a Cu 2+ 
ion in the top CuO layer of a current carrying HTS. This 
strong defect could either be a unitary impurity like Zn 2+ 
or simply a Cu-vacancy, and is well known to induce a 
near-zero-bias resonant peak (NZBRP) next to the site 
of the defect in a dSC without J s [3,4]. In the follow- 
ing, we investigate the LDOS images and their Fourier 
components for several values of the bias energy E and 
J s . In addition the LDOS at sites next to and far away 
from the impurity as a functions of E are calculated and 
their behaviors under various J s will be presented and 
discussed. 

The BCS Hamiltonian describing the impurity effects 
in a dSC carrying a supercurrent can be written as 
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where e k is the band structure of the d-wave supercon- 
ductor, p is the chemical potential to be determined by 
doping, V s is the on site potential of the nonmagnetic 
impurity located at the center of lattice, q s — (m * /2)v s 
with v s the supercurrent velocity, and m* the mass of a 
Cooper pair, A qs (k) = A qs cos(26l) is the superconduct- 
ing order parameter in the presence of J s , 9 is the angle 
between the wave vector k and the antinodal direction 
of the d-wave superconductor, and A qs is determined by 
the gap equation [7] 

Trln^ = J d0co S 2 (20)ln[ff (0) + ^g 2 (4>) - 1] , (2) 



(1) kp and Ep are the Fermi momentum and energy, re- 
spectively, <f> is the angle between q s and the antinodal 
direction, and the integraton in Eq. (2) is from to 2ir 
with the constraint g 2 (4>) — 1 > 0. The solutions of Eq. 
(2) with = and \ are presented in Fig. 1(a). In 
Ref. [7] , we also derived the thermodynamic critical cur- 
rents isc(0) = 0.238en«FAo at q = q c {0) = 0.35A and 
j sc (|) = Q.22henv F A at q = g c (f ) = 0.39A for su- 
percurrent J s along the antinodal and nodal directions, 
respectively [see Fig. 1(b)]. Similar results on the or- 
der parameter and the critical currents have also been 
obtained in Ref. [8]. 
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The Hamiltonian (1) is exactly soluble by the Bogoli- 
ubov transformation and the Green's function technique 
[10]. After a tedious but straightforward calculation, we 
obtain the expression for LDOS near a strong impurity 
(3) in the unitary limit (i.e. V s — > oo) 
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where N is the site number of lattice, and 
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iw n - 5(£k+q 3 - e-k+qj - (-l) y £ , qs (k)' G kl/ (q s ,iu;„) in Eq. (5)], which leads to different gaps 

for different momentum directions of a quasi-particle. 
This strongly modifies the LDOS and its Fourier 
components patterns. We shall see that the varia- 
tion of these patterns is sensitive to the supercurrent 



Obviously, when a supercurrent is appplied, the 
quasiparticle energy has a momentum-dependent shift 
|(ek+q s — £-k+qJ [sec the bare Green's function 
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applied and the bias voltage. However, they always 
have a reflection symmetry at an arbitrary energy if 
a supercurrent is applied along the antinodal or nodal 
directions. From Eqs.(4) and (5), we calculate the 
LDOS at several different energies and supercurrent 
velocities with the defect located at the center of a 
TV = 400 x 400 lattice. Here we adopt the band 
structure of Bi 2 Sr2CaCu208+,5 given by Norman et al. 
£k = — 0.5951 (cosk x + cosk y )/2 + O.lGSGcosfc^cosfcy — 
0.0519(cos2fcz + cos2fc y )/2 - 0.1117(cos2fc :c cos£;, / + 
cosk x cos2k y )/2 + 0.0510cos2fc x cos2fc J/ — /j, (eV) [11] 
, which corresponds to that of free electrons with 
k F = 1.639 and E F = 0.4203 cV, for the chemical 
potential /j, = —0.1238 eV for optimal doping (15%). 
Choosing A = 44 meV, and A q can be obtained from 
Fig. 1(a). 

Fig. 2 shows the 20 x 20 images of the LDOS p(r, w) 
at different energies and supercurrents with the impurity 
at its center. In order to understand the image patterns 
in Fig. 2, we plot the schematic Fermi surface of an op- 
timally doped HTS in the first Brillouin zone as shown 
in Fig. 3. In an STM experiment and when a quasi- 
particle is created near the Fermi surface at point O, this 
quasiparticle may be scattered elastically by the defect to 
other equivalent points (such as A, B, C, D, E, F and G) 
near the Fermi surface [12-16]. The wavevectors connect- 
ing O and the other points are referred as the modulation 
wavevectors and they are labeled as qA, <1b, etc., up to 
qp. If the point O is at the middle of the Fermi curve in 
Fig. 3, then the quasiparticle is at the nodal point. If O 
is moved to the zone boundary, then the quasiparticle is 
at the antinodal point. Because of the d-wave nature of 
the superconductivity, little energy is required to create 
a quasiparticle at the nodal point. But to create a quasi- 
particle at the antinodal point, a large bias energy in the 
order of the superconductivity gap is needed in an STM 
experiment. 

From Fig. 2, the LDOS at the impurity site vanishes 
regardless of the bias energies and the strength of J s , 
and it has the strongest intensity at the sites next to the 
defect when the quasi-particle energy (or the bias voltage 
times e) u = meV [see Fig. 2c(0)]. Near the impurity, 
the LDOS has a pattern of 4-fold symmetry with energy- 
dependent modulations in the absence of a supercurrent 
[Fig. 2a(0) to 2e(0)]. 

When \u>\ = and 16 meV, the resonant peaks still 
show up at (0, ±1) and (±1,0) and the modulation with 
the periodicity ~ 2a is along 45° from the Cu-0 bonds. 
For u> = meV, the point O is at the nodal point, and 
the modulation in Fig. 2c(0) comes from the wavevector 
qo = qE in Fig. 3. For \u\ = 16 meV, the pattern 
seems to be a result of combined contributions from qp 
and other modulation wavevectors along the directions 

(±7T,±7T). 

When |w| = 44 meV, The point O moves to the 
zone boundary or the antinodal point. The modulation 



wavevector qp = 2tt. This would give rise to x- and y- 
oricntcd (or Cu-0 bond oriented)) stripe-like structure 
with the periodicity ~ a in the LDOS around the im- 
purity, and this can be seen clearly in Fig. 2a(0). The 
LDOS here also exhibits a checkerboard pattern close to 
the impurity site due to the combined effect of the x- and 
y-oriented stripes. 

When a supercurrent J s is applid along the antinodal 
direction (i.e. = along the x-direction), the intensities 
of the resonance peaks on points (0, ±1) are higher than 
those on the points (±1,0) at \cu\ = and 16 meV. Near 
the critical current j sc (0), the LDOS developes a modu- 
lation perpendicular to the direction of supercurrent at 
\uj\ = 16 meV. When \u\ = 44 meV, the intensity of 
the modulation parallel to J s becomes smaller than that 
perpendicular to J s . When J s is applied along the nodal 
direction (i.e <f) — 7r /4 from the x-axis), the LDOS pat- 
terns only have minor changes except that some brighter 
spots near the defect site appear at \u>\ = 16 meV. We 
note that with increasing J Sl the maxima of the LDOS at 
|w| = and 44 meV are suppressed while those at \u\ = 
16 meV are enhanced. 

In order to further understand the supercurrent effects, 
we also calculate the images for the Fourier component 
of the LDOS (FCLDOS)( see Fig. 4). It can be easily 
seen that the influence of the applied suppercurent on 
the FCLDOS is more dramatic than on the LDOS image. 
Here some modulation wavevectors corresponding to the 
clastic scattering of quasiparticles from one point of the 
Fermi surface to another point as shown in Fig. 3 can be 
clearly identified in Fig. 4. 

When lo = meV and q s = 0, the modulation wavevec- 
tors qA = qc = If and qp, = qp due to the nodal quasi- 
particle scattering are clearly seen in the FCLDOS pat- 
terns [Fig. 4c(0)]. We note that the dip at qp.p has 
a strong intensity, which causes the LDOS to have a 
modulation along 45° to Cu-0 bonds [Fig. 4c(0)]. At 
the critical currents J sc (0) and J sc (j), the dips are su- 
pressed, but their positions seem not to shift [Fig. 4c(l) 
and 4c(2)]. 

However, at higher energy, the case becomes more com- 
plicated. At \ui\ = 16 meV, the dips at qA and qp arc 
clearly visible in the absence of supercurrent [Fig. 4b(0) 
and 4d(0)]. At four corners of the first Brillouin zone, 
there are four arcs due to the scatterings of the quasipar- 
ticles by the defect from one equal-energy banana countor 
(e.g. arcs OB in Fig. 2) to the opposite contour (e.g. arc 
DE in Fig. 3). It is these arcs that mainly produce the 
charge modulation along 45° from the Cu-0 bonds [Fig. 
2b(0) and 2d(0)]. At w = -16 meV, wc note that the 
peaks associated with qs are absent. Instead, four new 
dips at q2 show up. When a J s is applied, these peaks 
and arcs for |u>|=16 meV are supressed or enhanced, and 
even vanish near the critical current, but their position 
have little shift [see Fig. 4b(l), 4b(2), 4d(l) and 4d(2)]. 

When |w| = 44 meV, Fig. 4a(0) and 4e(0) show the 
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images of the FCLDOS associated with the scattering of 
antinodal quasiparticles. The peaks correspond to qa and 
qi are clearly seen here. The modulation wavevectors q\ 
are due to the superposition of those peak arcs induced by 
the scatterings of quasiparticles from one antinodal part 
of banana contour to the neighboring part. Obviously, 
the equal-energy banana contour becomes wide for J s 
along the antinodal direction while it shrinks to a point 
for J s along the nodal direction [Figs. 4a(l), 4a(2), 4c(l) 
and 4e(2)]. From Fig. 2 and Fig. 4, it can be seen clearly 
that the effect of J s on the FCLDOS is much pronounced 
than that on LDOS. 

In order to examine the change of modulation wavevec- 
tors with J s , we present the FCLDOS along the antin- 
odal and nodal directions at u = mcV and for several 
supercurrent strengths in Fig. 5. With increasing the su- 
percurrent velocity, the dips associated with qA,c,F are 
suppressed and finally disappear [see Fig. 5(a) and 5(c)]. 
However, the dips corresponding to qD,E are only sup- 
pressed for J s along the antinodal direction or (f> = 
while they are first enhanced, then weakened, and have 
a tiny shift for the J s along the nodal direction or <fi — f 
[Fig. 5(b) and 5(d)]. Similar results hold at higher en- 
ergy. We note that the dips at \q\ = 1.6ir in Fig. 5(a) 
and 5(c) cannot be induced by quasiparticle scatterings, 
which are also suppressed with increasing J s . We think 
that these dips at qA,c,F and qu,E are due to the mani- 
festation of quasiparticle destructive interference due to 
the sign change of the d-wave gap function on the Fermi 
surface. 

We have obtained the LDOS and FCLDOS induced 
by a strong defect such as a Zn impurity. Now we turn 
our attention to the STM experiments. In the STM ex- 
periments [5], a zero bias resonant peak was observed 
at the Zn sites. However, theoretical calculations give a 
vanishing LDOS at the impurity sites, contrary to the 
experimental observation (see Fig. 2). Because of a Bi 
atom in the top (BiO) layer and, more importantly, an 
O atom in the second (SrO) layer block the tunneling 
current coming from the STM tip to directly probe the 
impurity site [17], the experimentally observed LDOS at 
the impurity or Cu site should be approximately equal to 
the sum of those on four nearest neighbor sites around it 
[18], i.e. 

Pcxpt (R, u) p(R + S, uj) , (6) 

where 5 denote the nearest neighbor sites of the impurity 
or Cu ions. 

Taking into account this blocking effect, we present 
Pcxpt (R, w) curves at the impurity site (0,0) and the 
points (0,-1), (1,-1) and (2,0) for several J s along 
4> = and = 7r/4, respectively in Fig. 6. Obviously, 
the NZBRP on the impurity site and its neighbor sites 
are strongly suppressed and only slightly broadened with 
increasing J s in both directions. No splitting is clearly 



visible. On the other hand, the superconducting coher- 
ence peaks shows some suppression and splitting near the 
critical current along the antinodal direction with <j> = 
while their separation widens with increasing J s along 
the nodal direction with <fi = 7r/4. We further notice 
that the suppression of the NZBRP is insensitive to the 
direction of J s . A relevant work [20] studied the NZBRP 
due to a unitary impurity in the presence of a magnetic 
field and away from the vortex cores. The magnetic field 
effect was considered by including the Doppler shifts [21] 
in the energies of the quasiparticles and it thus generates 
circulating supercurrent in the sample which is similar 
but not identical to the case we studied. The NZBRP 
displayed in Fig. 1 of Ref. [20] and Fig. 6 in the present 
paper are both suppressed by the magnetic field B or J s , 
but it appears that their LDOS due to the impurity seem 
to lose a lot of spectral weight while ours practically re- 
mains as a constant as B or J s increases. Further work 
is needed in order to understand the difference between 
these two works. 

In summary, we have investigated the supercurrent 
effects on the impurity resonance states in d-wave su- 
perconductors. The LDOS and FCLDOS patterns in- 
duced by a strong impurity have a reflection symme- 
try if a supercurrent is applied along the antinodal or 
nodal directions. The suppression and broadening of the 
resonant peak and the superconducting coherence peaks 
are due to the anisotropic gap induced by a supercur- 
rent. Future STM experiments need to be performed 
in order to test these predictions. On the other hand, 
the midgap-surface-state-induced ZBCP in the tunneling 
conductance characteristics between a normal metal and 
a d-wave superconductor (dSC)[7] looks similar to the 
strong impurity induced NZBRP. But their dependences 
on the applied supercurrent J s are quite different. 

When strong defects such as Cu vacancies and micro- 
crystals with edges exposed to (110) direction are both 
present on the surface of a HTS sample, STM experi- 
ments should be able to distinguish them by analyzing 
the J s dependences of the zero-bias conductance peak 
induced by an isolated defect and that induced by the 
surface midgap states. 

The authors wish to thank Prof. S. H. Pan for help- 
ful Discussions. This work was supported by the Texas 
Center for Superconductivity and Advanced Materials at 
the University of Houston and by the Robert A. Welch 
Foundation (Ting). 



[1] C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 

(2000) and the references therein. 
[2] C. -R. Hu, Phys. Rev. Lett. 72, 1526 (1994). 
[3] A. V. Balatsky, M. I. Salkola, and A. Rosengren, Phys. 

Rev. B 51, 15547 (1995). 



5 



[4] M. I. Salkola, A. V. Balatsky and D. J. Scalapino, Phys. 

Rev. Lett. 77, 1841 (1996). 
[5] S. H. Pan, E. W. Hudson, K. M. Lang, H. Eisaki, S. 

Uchida, and J. C. Davis, Nature, (London) 403, 746 

(2000). 

[6] A. M. Zagoskin, Quantum Theory of Many-Body Sys- 
tems: Techniques and Applications (Springer, 1998). 

[7] Degang Zhang, C. S. Ting, and C.-R. Hu, Phys. Rev. B 
70, XXXXX (2004); cond-mat/0312545 

[8] I. Khavkine, H.-K. Kee, and K. Maki, |cond-mat /0405236 

[9] J. Ngai, Y. C. Tseng, P. Morales, V. Pribiag, J. Y. T. 
Wei, F. Chen, and D. D. Perovic, Appl. Phys. Lett. 84, 
1907 (2004). 

[10] C. Pepin and P. A. Lee, Phys. Rev. B 63, 054502 (2001). 
[11] M. R. Norman, M. Randeria, H. Ding, and J. C. Cam- 

puzano, Phys. Rev. B 52, 615 (1994). 
[12] J. E. Hoffman, K. McElroy, D.-H. Lee, K. M. Lang, H. 

Eisaki, S. Uchida, and J. C. Davis, Science 297, 1148 

(2002). 

[13] K. McElroy, R. W. Simmonds, J. E. Hoffman, D.-H. Lee, 
J. Orenstein, H. Eisaki, S. Uchida, and J. C. Davis, Na- 
ture (London) 422, 592 (2003). 

[14] Qiang-Hua Wang and D.-H. Lee, Phys. Rev. B 67, 
020511(R) (2003). 

[15] Degang Zhang and C. S. Ting, Phys. Rev. B 67, 
100506(R) (2003); Phys. Rev. B 69, 012501 (2004). 

[16] Lingyin Zhu, W. A. Atkinson, and P. J. Hirschfeld, Phys. 
Rev. B 69, 060503(R) (2004). 

[17] Q. Wang and C.-R. Hu, unpublished. 

[18] J.-X. Zhu, C. S. Ting, and C.-R. Hu, Phys. Rev. B 62, 
6027 (2000). For an alternative model, see Ref. 19. 

[19] I. Martin, A. V. Balatsky, and J. Zaanen, Phys. Rev. 
Lett. 88, 097003 (2002). 



[20] K. V. Samokhin, Phys. Rev. B 68, 104509 (2003). 
[21] G. E. Volovik, Pis'ma Zh. Eksp. Theor. Fiz. 58, 457 
(1993). [JETP Lett. 68, 469 (1993)]. 

Fig. 1: Dependences of the superconducting order pa- 
rameter on the normalized supercurrent-velocity param- 
eter q for a d-wave superconductor (a) and the corre- 
sponding dependences of supercurrent density on q (b). 

Fig. 2: The 20 x 20 images of the LDOS p(r,u>) at 
uj = —44 meV, —16 meV, meV, 16 meV and 44 meV 
(from top to bottom) and q = 0,q c (0) and q c (j ) (from 
left to right) for a unitary impurity at its center. 

Fig. 3: Schematic Fermi surface of high-Tc cuprate 
superconductor. The modulation wave vectors connect- 
ing different points of the Fermi surface with the same 
energy gap are shown in the absence of a supercurrent. 

Fig. 4: The FCLDOS at u = -44 meV, -16 meV, 
meV, 16 meV and 44 meV (from top to bottom) and 
q = 0, Qc(0) and q c (j) (from left to right) in the first 
Brillouin zone for a unitary defect. 

Fig. 5: The FCLDOS along the antinodal and nodal 
directions at u> = meV and different supercurrents for 
a unitary defect. 

Fig. 6: The predicted, blocking-model-corrected 
LDOS p oxp t(R, to) at the sites (0, 0), (0, -1), (-1, -1) and 
(2, 0) (from top to bottom) for supercurrents along the 
antinodal (left) and nodal (right) directions, when a uni- 
tary impurity is located at the (0,0) site. Solid: q = 0, 
dash: q = 0.2A and dot: q = q c (4>)- 
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